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C^ ■ Abstract 

The standard Regge poles are considered as inputs for two unitarization 
C^^ '_ methods: eikonal and U-matrix. It is shown that only models with three in- 

put pomerons and two input odderons can describe the high energy data on pp 
r^ , and pp elastic scattering including the new data from Tevatron and LHC. How- 

ever, it seems that the both considered models (eikonal and [/-matrix) require a 
Oh' further modification (e.g., to explore nonlinear reggeon trajectories and/or non- 

(— I . exponential vertex functions) for a more satisfactory description of the data at 

19.0 GeV< ^ < 7 TeV and 0.01 < |t| <14.2 GeV^. 

(N 
> 
O 1 Introduction 

00 . 

O ■ The first results of the TOTEM experiment on elastic proton-proton scattering at 7 

■^ I TeV were published in 2011 [T]. Surprisingly none of models (most successful at lower 

^ I energies and wide intervals of t) could correctly predict the t-dependence of the pp 

differential cross sections at LHC energy [21 13]. 

The data on elastic pp and pp scattering at energies up to Tevatron one were quite 

successfully described in three types of models. The fist type deal with the n intercept 

larger than one pomerons (e.g. Donnachie-Landshoff model, see [5] and the recent 

papers [6], [?]). However such models have defect: they violate the Froissart bound 

on the total cross sections and have to be unitarized even if they are in agreement 

with the data. This fearture give rises to various methods of unitarization of the 

input amplitudes violating the unitarity. The most known models of such type are 

the eikonal (see a general discussion on the eikonal method in [S]) (or quasi-eikonal 

[H]) and [/-matrix [lOJ (or its generalization [H]). The models of the third type are 

constructed taking into account restrictions of amplitude analyticity and unitarity from 

the beginning. 
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The TOTEM data have been analyzed in a few recent pubhcations. Here, for the 
sake of shortness we do not review the obtained results. The corresponding references 
and discussion on various models and its agreement or disagreement with the old 
data and recent DO [1] and TOTEM [2] ones are well presented in [12] and in the 
comprehensive review [13] . 

We will concentrated here on the eikonal and [/-matrix models in their "standard" 
form with the simplest input reggeon amplitude. The name "standard" refers to em- 
ploying reggeons with linear trajectories and exponential form of vertex functions. Our 
aim is to check if such models are able to describe data on pp and pp elastic scattering 
at high energies and in all available transferred momenta including the new data . We 
will also to determine the minimal number of the input reggeons that contribute to 
amplitude and provide a reasonable agreement with the data. 

2 "Standard" schemes of unitarization 

2.1 Impact amplitudes 

Let define the elastic amplitude in impact parameter (two-dimensional vector b) rep- 
resentation, or in short impact amplitude, as 

H{s, h) = ^J ^e^%(., t) = ^J dkkMb^/^t)A{s, t) (1) 



with the inverse transformation 

_, oo 

A{s, t) = levr^s f -|-^e-*%(s, b) = 8ns f dbbJoibV^)H{s, b) (2) 



where A[s, t) is an elastic amplitude and Jo{z) is the Bessel function, k is two-dimensional 
vector and t ~ — fc^. 

One can obtain from the unitarity equation for A{s,t) (see e.g. [^) the unitarity 
equation for H{s, b) which in a general form can be written as 

'^smHis^b) = \H{s,b)\^ + G,nei{s,b), (3) 

where Ginei{s,b) > takes into account inelastic intermediate states from the original 
unitarity condition for A[s,t). Therefore 

3mi7(s, b) > 0. (4) 

and from the unitarity one can obtain 

\H{s,b)\<l. (5) 

For the observed cross sections one can derive the following expressions 

oo 

at{s) = -SmA(s, 0) = Svr / dbb'^mH{s, b) (6) 
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cTei{s) = j^ J dt\A{s,t)\' = 8n J dbb\H{s,b)\^ (7) 



-oo 

oo 



cTineiis) = SvT f dbb{^mH{s,b) - \H{s,b)\''). (8) 







2.2 Output amplitudes 

We consider two known methods for unitarisation of input amplitudes h°-^{s,b) where 
a = p OT a = p. Both of them can be considered as the particular cases of sum of 
multireggeon exchanges (see fig. [1]). Just for illustration we show here only pomeron 
contributions, a more realistic case will be considered later on. 





Figure 1: Multipomeron contributions to elastic scattering amplitude 






Figure 2: Amplitude of interaction of two hadrons with n pomerons in a pole approx- 
imation but with phenomenological factor Xa{n). 



Assuming that the two-hadrons-n-pomeron amplitude is proportional to the prod- 
uct of the two-hadron-pomeron vertices depending only on n (it is a pole approximation 
for intermediate states) as shown in fig. HI one can obtain 

n=l 

Moreover, assuming either A(n) = A" or \{n) = A"n!, where A = AaAp(= Xap), 
we obtain from Eq. i^ two well-known schemes of pomeron unitarization: eikonal 
[SI [HI [15] or quasi-eikonal [S] and [/-matrix or quasi- [/-matrix models [101 ITT| . 
Eikonal and Quasi-eikonal unitarisation 

2iXaph''P{s,b) _ 1 

^"'^''^^ = 2^A ' VH = (Vr (10) 

If Xap = 1 then Eq. (llOj) describes the pure eikonal unitarization. 



U-matrix and quasi- U -matrix unitarization 

For \ap = 1/2 we have the pure U-matrix unitarization. 
In the both cases 

Aap > 1/2 (12) 

because full amplitude must satisfy the inequality \H°'^\ < 1. 

2.3 Standard Regge pole input amplitudes 

For input amplitudes we consider the standard Regge pole contributions 

<(^'^) = E«+(^'^)±E^-(^'^) (13) 

R+ R- 

At high energy {s ^ mV) the input amplitudes can be written as following 

a±(^,t)=("/)^iJ±W(-^Sr±^^ (14) 

with s = s/sq, so=lGeV^. 

There are two kinds of Regge poles contributing to amplitudes. The first ones are 
poles with intercepts close to 1 (or large than 1). They are crossing even pomeron (or 
pomerons) and crossing odd odderon (or odderons). The second kind of contributions 
are so called secondary reggeons with intercepts less than 1. They are f,a2,uj,p with 
a{0) ~ 0.5 and other with lower intercepts. For high energy pp{pp) scattering it is 
sufficient to consider in the amplitudes one effective crossing even secondary reggeon 
(-R+) and one effective crossing odd secondary reggeon (-R_). 

In what following we keep the arguments of the [T^ and consider the secondary 
crossing even and crossing odd reggeons with the fixed intercepts and slopes (strictly 
speaking, effective trajectories are not coincided with those for /- and u- reggeons 
while at high energy their contributions to amplitudes are closed) 



a+(0) = 0.69, aV = 0.84GeV-2 
a_(0) = 0.47, a'_ = 0.93GeV 



/ _ n nQn„A/-2 (15) 



3 One (two) pomeron and one odderon model 
(lP(2P)+0-model) 

A. One exponent form of the vertex functions gpi{t),i = 1,2 and go{t)- 

For eikonal form of the amplitude fITT]) one and two pomeron models with the 
exponential vertex functions {gi{t) = giexp{Bit)) input amplitudes have the form 

a,{s,t) = 77,(7, exp(5,t)(-^§)"'(0)+"^* (16) 



where i = Pi,P2,0,+,— and 7]i = —1 for pomerons and crossing even secondary 
reggeon, rj = i for odderon and crossing odd reggeon. In accordance with Eq. ([T]) the 
impact amphtudes hi{s,b) are the following 

h,{s, h) = ry,^^(-^S)°^(°) exp(-6V4r2) (17) 

where 

r2 = 5i + a;in(-z5), z = Pi, P2, 0, +, -. (18) 

These models have been considered by Petrov and Prokudin [l4j. They found that the 
P+0 model as well as 2P+0 one, the both with linear trajectories and with exponential 
form of the vertex functions g{t) (for all reggeons including pomerons and odderon), 
failed to describe data with a sufficient quality. 

B. Two exponents form of the vertex functions gpi(t), i = 1,2 and go{i)- To check 
if a more complicated vertex function g{t) can improve the data description we have 
considered for pomeron and odderon terms (or for some of them) a combination of two 
exponents with the different slopes B. Namely, we write for i = Pi, P2, O 

9i{t) = gi[ciexp{Biit) + (1 - Ci) exp{Bi2t)]. (19) 

Changes for the impact amplitudes hpi{s,b) and ho{s,b) are quite evidently follow 
from the Eqs. flTTini)]) . 

We do not give here the results of the fit for this modified model, we rather note 
that such an extension of the model in spite of some its improvement does not lead to 
a satisfactory description of the TOTEM data. 

4 Three pomeron and one odderon model 
(3P+0-model) 

This model has been suggested and considered by V. Petrov and A. Prokudin in 
|14j . The input amplitudes h°'''{s,b) in the model are written as contribution of three 
pomerons, one odderon and two secondary (even and odd) reggeons 

h%{s,b) =hp,{s,b) + hp^{s,b) + hp^{s,b) + h+{s,b) , . 

±{hois,b) + h4s,b)) ^^^^ 

where input impact amplitudes are given by Eq. flTTl) . Then the eikonalized amplitudes 
APP{s,t) are defined by Eq. fITOl) with A"^ = 1 and Eq. ([2]). They have obtained very 
good description of the pp and pp data for the total cross sections crtot{s), the ratios of 
real to imaginary part of the forward elastic scattering amplitudes p{s) at ^/s > 8 GeV 
and for the differential cross sections da{s,t)/dt at ^/s > 23 GeV and 0.01 < |t| < 14 
GeV^ (this data set contains totally 961 points). 

It is necessary to note here that our data set differs of those used in [TH |T5] . Our 
data set has been proposed by J.R. Cudell, A. Lengyel and E. Martynov (CLM) [16], 
who built a coherent set of all existing data for 4 < ^/s < 1800 GeV and < |t| < 15 
GeV^ JW\ which can be found in the HEP DATA system [T71 |T2] . CLM have performed 



a detailed study of the systematic errors and gathered into a common format more than 
260 subsets of data from more than 80 experimental papers. We suggest to use the CLM 
data set as a standard dataset. The latest (with some corrections) updated version of 
it including TOTEM [H El |3] and DO [1] data is available online [19]. The set used for 
the given analysis contains 1723 points In the region described above. There are only 
31 points from the 3 groups, 8 points at a/s =26.946 GeV, 11 points at 30.7 GeV and 
12 points at 53.018 GeV for dapp/dt were excluded from the all data presented in |TB] 
because these groups are strongly deviated from the rest data points and can shghtly 
distort the fit. 

As was shown in the paper [1] where the first data at LHC energy 7 TeV being 
presented, the model [H] does not describe the TOTEM data. This statement is 
confirmed in the recent paper [IS] -by authors of the model (we reproduce Fig. 4 of 
their paper in the fig. [5a|) . We have performed refit of the model with our data set 
including the TOTEM data and found that the resulting theoretical curve for da/dt 
at 7 TeV is still out of the data at |t| > 0.4 GeV^ as shown on the fig. [3bl Thus we 
can conclude that the 3P+0-model should be modified either adding new terms or 
using nonlinear pomeron (odderon) trajectories and more complicated non exponential 
vertex functions. In what follows we consider the first possibility, namely we would 
like to check if an adding the second odderon term can help to improve agreement of 
the model with new data. 



E 



\ '^ 


j,^|%l|f||W^^ 




\ 


"(GeV^)' 


\ pp TOTEM 


\ ^"""^^ 







-t (GeV^ 




1.5 2 

|t| (Gev2) 



(a) from the 



(b) after refitting 



Figure 3: Description of the TOTEM data in the three pomerons and one odderon 
model [H] 



5 Three pomerons and two odderons model 
(3P+20-model) 

In this section we consider the minimal extention of the 3P+0-model ( 120|) adding the 
second odderon term. 



h%{s,h)=hp{sM + K{sM±{ho{sM + h-{s,h)), 

3 2 

hp{s,h) = T.hpXs,b), ho{s,b) = E^o,(s,6) 



(21) 



j=i 



i=l 



where each input term is still in a "standard" form (ITi|) . (IT6|l with linear trajectory 
ai = ai(0) + a^t and exponential vertices gi(t) = giexp{Bit). 

Later we will check another possibilities to improve the eikonal and U-matrix models 
(nonlinear trajectories, non exponential vertices and so on). 

We considered the both methods of unitarization, eikonal one, Eq. (llOj) . with A'^^ = 1 
and [/-matrix one, Eq. flTTl) . with A'*^ = 1/2. Fit has been performed with the data in 
the region 

for crtot(s) and p(s) at -\/s > 5GeV, 
for da(s,t)/dt at ^i > 19.0GeV (22) 

0.01 < |t| < 14.2GeV^ 

The resulting description of the data is shown in the fig. HJ-fig. El Parameters of the 
models are given in the Table [TJ 




(a) Total pp, pp cross sections. 



150 


- 


pp cros 


s sections 


140 


- 


X 


pp data 
pp data 


130 


- 





TOTEM 


120 


- 


n 
o 


ALICE 
ATLAS 


110 


- 


i> 


CMS 


100 


- 




— Eikonal 




Eikonal 









- U-matrix 


BO 


- 





- U-matrix 


70 


- 





- U-matrix 




•Js (GeV) 

(b) Total, elastic and inelastic pp cross 
sections. 



Figure 4: Total, elastic and inelastic cross sections in 3P-f20 models. Solid (dashed) 
line shows the eikonal ([/-matrix) results. 

Let us briefly comment the results. 

Addition of the second odderon term leads to certain improvement of the data 
description. 

Agreement with the TOTEM data is essentially better but some defects are ap- 
peared at lower energies, especially at large \t\ and for pp differential cross sections (see 
figs. [6] and [7]). 
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Figure 5: The ratio of the real to imaginary part of the forward scattering amphtude. 
Sohd (dashed) hne shows the eikonal ([/-matrix) model results. 
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Figure 6: Elastic pp and pp scattering in the 3P-I-20 model. Solid (dashed) line shows 
the eikonal (^/-matrix) results. 



Big difference for ppp and ppp (shown in fig. [5]) is disappeared at very high energy. 
In the eikonal model the difference is negligible at ^/s > 10'' GeV while in the U-matrix 
model it is at a/s > 10^ GeV. Thus an asymptotic regime comes quite late in energy. 

In the Table [2] the calculated values of the total, elastic and inelastic cross section 
are given for the LHC energies. One can see a very good agreement with the results of 
TOTEM [3], ALICE |22], ATLAS [20] and CMS |2l]. 
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Figure 7: Description of the TOTEM data and predictions for high |t| of the 3P+20 
models. Solid (dashed) line shows the eikonal ([/-matrix) results. 

More problems are for the U-matrix model rather for the eikonal one, however both 
models do not strongly contradict to the TOTEM data. 

6 Conclusion 

Answering on the question put in the title we would like to say that in our opinion 
the obtained results demonstrate an ability of the unitarized eikonal and U-matrix 
3P+20 models to describe the data up to LHC energies, however, it seems that for 
further improvement nonlinear and/or non exponential (for vertices) effects are very 
important. Exploring another possibility, namely, construction of the model that does 
not explicitly violate unitarity (at least the unitarity bounds on the cross sections) just 
from the beginning its construction is in a progress, results will be presented later. 
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Table 1: Parameters of the eikonal and [/-matrix models obtained by fitting to the 
data. Parameters a^ and Bi are given in GeV~^, other parameters are dimensionless. 
Errors are taken from the MINUIT output, 'x^ /dof = x^/i^p — Npar) where Np is the 
number of experimental points in the fit, Npar is the number of free parameters in the 
model. 



Parameters 


Eikonal model 


U-matrix model 


value 


error 


value 


error 


api 


1.2209 


0.0012 


1.2654 


0.0029 


< 


0.1171 


0.0007 


0.1041 


0.0005 


9Pi 


1.3249 


0.0196 


0.1371 


0.0037 


Bp, 


0.3977 


0.0055 


0.0 


0.001 


ap2 


1.1675 


0.0006 


1.1208 


0.0006 


«k 


0.3036 


0.0014 


0.2324 


0.0048 


9P2 


8.753 


0.064 


26.879 


0.142 


Bp2 


0.4125 


0.0089 


5.325 


0.0455 


aPs 


1.0703 


0.0005 


1.089 


0.001 


«k 


0.5912 


0.0024 


0.2393 


0.0017 


SP3 


49.138 


0.198 


23.657 


0.161 


Bp, 


1.2046 


0.0169 


1.2210 


0.0114 


aoi 


1.220 


0.0004 


1.1889 


0.0001 


»'o. 


0.0723 


0.7E-05 


0.0422 


0.00004 


90i 


-27.133 


0.007 


-35.072 


0.026 


Bo, 


0.01384 


0.00003 


1.052 


0.0003 


ao2 


1.2196 


0.00002 


1.185 


0.0001 


»'o. 


0.07222 


0.5E-05 


0.04197 


0.00004 


902 


27.171 


0.0067 


36.969 


0.027 


B02 


0.01374 


0.00003 


1.0671 


0.0003 


a+ 


0.69 


fixed 


0.69 


fixed 


«v 


0.84 


fixed 


0.84 


fixed 


9+ 


221.55 


1.18 


244.49 


1.05 


B+ 


5.460 


0.175 


0.0 


0.004 


a_ 


0.47 


fixed 


0.47 


fixed 


a'_ 


0.93 


fixed 


0.9 


fixed 


9- 


149.17 


4.08 


182.98 


2.85 


B 


7.057 


1.086 


1.337 


0.189 


xVdof 


1.673 


1.973 



Table 2: Cross sections and ratio of real to imaginary part of the forward scattering 
pp amplitude in the eikonal model 3P+20 defined by the Eqs. flTOj) . fl2T]) . 



V^{TeV) 


atot (mb) 


aei (mb) 


ainei (mb) 


P 


7 


98.97 


25.53 


73.45 


0.153 


8 


101.5 


26.48 


75.02 


0.152 


13 


111.14 


30.18 


80.96 


0.150 


14 


112.68 


30.78 


81.9 


0.150 
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